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ABSTRACT 

We report on a confirmed galaxy cluster at z = 1.62. We discovered two concen- 
trations of galaxies at z 1.6 in the Subaru/XMM-Newton deep field based on deep 
multi-band photometric data. We made a near-IR spectroscopic follow-up observation 
of them and confirmed several massive galaxies at z = 1.62. One of the two is associated 
with an extended X-ray emission at 4.5(7 on a scale of 0'.5, which is typical of high-z 
clusters. The X-ray detection suggests that it is a gravitationally bound system. The 
other one shows a hint of an X-ray signal, but only at l.Scx, and we obtained only one 
secure redshift at z = 1.62. We are not yet sure if this is a collapsed system. The 
possible twins exhibit a clear red sequence aX, K < 22 and seem to host relatively few 
number of faint red galaxies. Massive red galaxies are likely old galaxies - they have 
colors consistent with the formation redshift of =3 and a spectral fit of the brightest 
confirmed member yields an age of l.^tol Gyr with a mass of 2.51q'^ x 10"*^^ Mq. Our 
results show that it is feasible to detect clusters at z > 1.5 in X-rays and also to per- 
form detailed analysis of galaxies in them with the existing near-IR facilities on large 
telescopes. 



Subject headings: galaxies: clusters: individual — galaxies: luminosity function, mass 
function — X-rays: galaxies: clusters 
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1. Introduction 

The galaxy evolution is closely linked to the structure formation of the Universe. Isolated 
galaxies and those in clusters evolve in different ways and the differential evolution results in 
the strong environmental dependence of galaxy properties observed in the local Universe. The 
environmental dependence of galaxy properties is already strong at z = 1. Intensive studies of 2; ~ 1 
clusters have shown that red galaxies have already be c ome the dominant population in clus ters at 



1 (e.g., iBlakeslee et al.ll2003l : iNakata et al. 



2005 



Lidman et al.ll2008l : iMei et al.ll2009l ^. The 

spectroscopically coii firmed highest redshift X-ray cluster known until now is located at z = 1.4 5 



(IStanford et al.l 120061 ). but red galaxies are abundant even in this highest-z cl uster (IHilton 



et al 



2OO9I ). A few authors reported a possible clusters at z > 1.5. For example, Kurk et al. ( 20091 ) 
presented an over-density of galaxies at z = 1.6, which may la ter collapse to a clust er. It is not 



detected in X-rays down to a limit of < 1 x 10~ ergs s~ cm~ . lAndreon et al.l (|2009l ) presented a 
cluster at Zphot ~ 1-9, but its redshift has been questioned by Bielby et al. (in prep) — it is likely 
a complex system with one confirmed at z = 1.1 and another one possibly at z ~ 1.5 (see Bielby et 
al. for details). It is challenging to identify clusters at z > 1.5, but one has to explore this redshift 
regime to identify the epoch when red galaxies become the dominant population in clusters and 
how the environmental dependence of galaxy properties is established. 

We are conducting a distan t X-ray cluster su r vey in the Subaru/XMM-Newton Deep Field 
(SXDF). We refer the reader to iFinoguenov et al.l (|2009bl ) for details of our survey and the con- 
struction of the deep multi-band photometric catalog. As part of the survey, we identified two 
concentrations of red galaxies first by their red sequence and then we found one of the sources is 
also securely detected in X-rays. We estimated their redshifts to be z ~ 1.6 both from the location 
of the cluster red sequence and from photometric redshifts. Not only they have very similar red- 
shifts, but they are also close to each other on the sky; they are separated only by ~ 2.5 arcmin 
(~ 1.3 Mpc at 2; = 1.6). Followed by the initial photometric identification, we carried out a near-IR 
spectroscopic foll ow-up observatio r i of th e possible twin clusters and we report on the results in the 
paper. Recently, iPapovich et al.l (j20ld ) presented a completely independent study of one of the 
clusters. Readers are referred to their paper for a similar, but independent analysis. 

The layout of the paper is as follows. We first describe the near-infrared spectroscopic follow- 
up observation and present a discovery of the most distant X-ray cluster from the observation in 
section 2. Section 3 discusses our results and summarizes the paper. Unless otherwise stated, we 
adopt Hq = 72km s~^ Mpc~^, Qm = 0.25, and Q\ = 0.75. Magnitudes are on the AB system. 



We note that we have revised the X-ray catalog presented in IFinoguenov et al.l (|2009bl ) by adding the 0.3-0.5 keV 
band data from XMM. We extract the images separately using the single events, produce a corresponding background 
estimation, and add them to the mosaic only at the very last stage. This increases signals from distant low-mass 
clusters. 
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2. Observation and Discovery 



We carried out near-IR multi-obje ct spectroscopy witii MOIRCS on the Subaru telescope 
(jlchikawa et al.ll2006l : ISuzuki et al.ll2008l ) on the 3rd-4th of December 2009. We selected targets for 
spectroscopy using photometric redshifts ; the first priority is given to red galaxies at 1.4 < Zpj^^ot < 
1.8, and second priority to blue galaxies at the same redshift range. The observing conditions were 
poor; 1 — 1.2 arcsec seeing in the J-band most of the time and we occasionally had clouds and high 
humidity. We used the z J500 grism with a slit width of 0.8 arcsec, which gave a resolving power of 
R ^ 500. The exposures were broken into 15 min. to 20 min. each and the total integration time 
amounted to 320 mi n. The data were reduced in a standard manner with custom-designed pipeline 
(jTanaka et al.ll2009l ). The spectra were visually inspected and redshifts and confidence flags were 
assigned to each object. We obtained 11 secure redshifts and 5 possible redshifts out of 39 observed 
objects. We present sample spectra in Fig. [TJ 

Fig. [2] summarizes our effort. There are two distinct concentrations of red galaxies. The 
concentration on the top-right (SXDF-XCLJ0218-0510) is detected in X-rays at 4.5ct, which suggests 
that it is a gravitationally bound system. The X-ray is extended towards West, but we suspect 
that this is because a gap between the pn detectors passes the Eastern part, where only MOS 
data contributed for the detection. It has two very bright galaxies and they have secure redshifts 
at z = 1.625 and z = 1.634, respectively. There are two more galaxies with secure redshifts at 
z ~ 1.625 in this system. Together with the clear concentration of the photo-z selected galaxies, 
we are highly confident that this system is a real cluster. This is the first confirmation of a distant 
cluster based on near-IR spectroscopy. 



Papovich et al.l (120101 ) confirmed several low-mass blue galaxies at the same redshift, o ne of 
which is a common object at z = 1.649. Including redshifts from iPapovich et al.l (|2010l ). we 



1.6230;°;°°°^ and a 



have 9 secure redshifts in this system, and we measure its redshift to be z 
velocity dispersion of 537I213 ^™ using the biweight estimator and gapper method, respectively 



(jBeers et al.lll990l ). We note that biweight estimator gives a velocity dispersion of 109^3q'^ km s 



being consistent with the one from gapper within the erro rs. We take the o ne from gapper in the 
following because of its robustness against poor statistics (IBeers et al.lll990l^. The virial mass and 



radius are M200 = l-ltoi x IO^^Mq and i?2oo = 0.54+^;^^ Mpc (ICarlberg et al.lll997l'). We note tha t 



+0.08 



one obtains a biased high velocity dispersion when only blue galaxies are used (IBiviano et al.ll2006l ). 
but half of the galaxies used here are massive red galaxies. Two galaxies at z = 1.642 and z = 1.649 
might be background galaxies because they increase the velocity dispersion to 1183^339 km s^^, 
which seems too large, but we are not sure at this point. 

The properties of the system derived from X-rays are consistent from those from optical. We 
detect an appar ent flux of 8.1 ifc 3.7 x 10 



-16 



ergs s 



-1 —2 

cm 



at 0.5-2.0 keV (37 counts) after the point 

1 



source removal ( Finoeuenov et al ] l2009bl ). which can be translated into Lx = 3.4=fc l.6x 10 ergs 
at 2: = 1.6 in the rest-frame 0.1-2.4 keV. Then, assuming the relation derived by iLeauthaud et al 
(j20ld ). we obtain a cluster mass of M200 = 5.7 ± 1.4 x IO^^Mq and a virial radius of -R200 = 0.440 
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Mpc, being consistent with the optical estimates. Clust er temperature is estima ted to be Tx 



1.7 lb 0.3 keV from the Lx — Tx relation as mentioned in lLeauthaud et al.l (|20ld ) 



We emphasize that most previous cluster studies did not remove fluxes from unresolved point 
sources (AGNs) once the extent of a cluster has been established. This procedure no longer works 
in deep surveys because the contamination from point sources is significant. We have therefore 
adopted the procedure in which we always remove the flux expected from detected point sources 
based on the shape of the PSF. In fact, we have attributed a large fraction (~ 50%) of the emission 
around the cluster to point sources. At the XMM exposure around the cluster, if a point source 
contributes substantially to the overall emission (> 30%), it can be detected individually. We 
can therefore exclude a possibility of large contamination from undetected point sources to the 
measured extended flux. However, some of the detected point sources might be associated with the 
cluster core and high resolution X-ray imaging will allow us to identify potential confusions. 

The one on the bottom-left (SXDF-XCLJ0218-0512) in Fig. [2] hosts an outstandingly bright 
galaxy, which is likely a cD galaxy of this system. We obtain a possible redshift of z = 1.625 for 
this galaxy. Interestingly, we observe a hint of an X-ray emission right on top of the cD galaxy 
at 1.5cT. Note that we do not observe nebular emission lines from the cD galaxy (see the bottom 
panel in Fig. [T]), but we need deeper observations to confirm if the X-ray signal is real and if it 
is from an AGN or the cluster core. We confirm another galaxy at 2; = 1.627 nearby. Although 
the concentration of galaxies around the possible cD galaxy is clear, we are not yet sure if this is a 
collapsed system due to the weakness of X-ray emission and to too few spectroscopic redshifts. For 
the rest of the paper, we assume that the two systems are real clusters at the same redshift and 
combine them to gain statistics. We note that our results remain unchanged if we use the secure 
cluster only. 



Results and Discussion 



The cluster red sequence is a ubiquitous feature of galaxy clusters. We draw color-magnitude 



diagrams in Fig. [3]using galaxies at < 1 Mpc from the twins. We refer the readers to lPapovich et al 



(j2010l ) for a similar analysis. The photo-z selected galaxies form a clear red sequence ai K < 22, 
while the sequence is not very clear at fainter magnitudes . We model the location of the red 



sequence with the updated version of the lBruzual &: CharlotI (j2003l ) code, w hich takes into a ccoun t 
eff'ects of thermally pulsating AGB s tars, usiiig the procedure described in iLidman et al.l (j2008l ). 
Here we assume the Chabrier IMF (jChabrieri |2003| ). A biweight fit to the red galaxies gives a 
sequence very close to the zj = 3 model sequence, suggesting that they are relatively old galaxies. 
The cosmic time between z = 3 and 1.6 is 2.0 Gyr. We note that we obtain Zf = 3 or higher if we 
use the z — J color as used in lPapovich et al.l (|20ld ) who obtained Zf = 2.4. We do not know the 
cause of the difference in at this point, but it could be due to the rather old UKIDSS catalog they 
used. Their paper is still under review at the time of this writing and we do not pursue this point 
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further. The z — K color is more sensitive to Zf at this redshiflj^ and we prefer to use it in this paper. 
As a sanity check, we perform a statistical subtracti on of fore-/background galaxies without using 
photo- 2: following the recipe by iTanaka et al.l (|2005l ) in the right panel. A clump of red galaxies 
is clearly visible in the right panel as well. Again, the sequence vanishes sX K > 22. There are 
a number of proto-clus ters at even higher redshi ft that have a confirmed over-density of low-mass 
star forming galaxies ()Miley fc: De Breuckl l2008l ) . but Fig [3] shows that a near-IR spectroscopic 
observation is essential to confirm the dominant population of massive galaxies in 2; > 1.5 clusters. 

To further quantify the red sequence, we plot luminosity functions (LFs) of red galaxies in Fig. 
m Here we define red galaxies as those having A|z — K\ < 0.2 from the best fitting red sequence. 
We confirmed that a small change in the definition does not alter the result below. But, we should 
note that the current statistics is very poor. We use galaxies at 1.4 < Zp^ot < 1-8 and statistically 
subtract the remaining contamination using the entire SXDF field as a control fi eld sa mple. We 
compare the z = 1.6 LF with those at lower redshifts taken from iTanaka et al.l (120081 ). We plot 
LFs relative to at each redshift to correct for the passive evolution. The LF of the red galaxies 
in the twins is similar to the group LF at z = 1.1, which suggests a deficit of faint red galaxies. 
On the other hand, it seems that clusters at z = 1.6 already host abundant massive galaxies. As 
seen in Fig. [3l these massive galaxies are mostly red galaxies and the clusters host few bright blue 
galaxies. It is hard to generalize the trend we see here with only a few clusters, but the implication 
would be that massive red galaxies might have become the dominant population in clusters already 
at z = 1.6. 

Finally, we present a spectral analysis of the brightest red galaxy with secure redshift in 
Fig, m We fit the observ ed spectra with model spectra generated by the updated version of the 
Bruzual &: CharlotI (j2003l ) code. We assume the exponentially decaying star formation histories 
and use age, dust extincti on, star format ion time scale, and metallicity as free parameters. We 
adopt the Chabrier IMF (jchabrierl bood j^ and fix the redshift at 

Zgpec- The fit gives an age of 



LS^Qo Gyr and a small extinction of Ay 



0.2^0:2 mag. 



Its star formation time scale is rather 



short, O.l^Q^ Gyr. Assuming the exponential decay, it assembled 80% of its stars at Zf ~ 2.5 that 
it would have at z = 0, which is in agreement with that inferred from the z — K color. Note that 
Zf = 2> obtained from the red sequence gives an average formation redshift of the red galaxies, while 
we are specifically discussing the brightest red galaxy with secure redshift here. The current star 
formation rate is < O.IMq yr~^, which is consistent with the absence of any emission lines in the 
observed spectrum. Its mass is estimated to be 2.51q'^ x IO^^Mq, revealing the presence of such 
a massive galaxy in a 2; = 1.6 cluste r. As seen in lower r edshift clusters, massive cluster galaxies 
form early and evolve passively (e.g.. iThomas et al.ll2005l ). 



This study is the first result on cluster galaxies at 2 > 1.5 using near-IR spectrograph. Our 



^ The z — K color of the model red sequence formed at zf = 5 is redder than the Zf = 2 sequence by 0.74, while 
the difference is only 0.26 for the z — J color. Thus, the z — K color is more sensitive to 2;/ than z — J. 

^ If we use the Salpeter IMF, we obtain 1.8 times higher stellar mass. The other properties are nearly the same. 
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observation suggests that z = 1.6 clusters exhibit a red sequence at bright magnitudes and are 
Ukely luminous in X-rays, which may have an implication for future high-z cluste r surveys. The 
cluster we discovered will evolve to a 4-5 keV cluster at z = (jvan den Boschll2002l ). Such clusters 
are a sensitive probe of cosmology at high redshifts and yet its estimated temperature preclude 
a detection in a Sunyaev-Zeldovich observation. This demonstrates the power of X-ray observa- 
tions in finding distant clusters. The redshift of z > 1.5 is the epoch when early dark energy is 
important and la rge-scale structures at z > 1.5 are a robu s t tracer of the primord ial local non- 
Gaussianity (e.g.. iBartelmann et al.ll2006l : ISadeh et al.ll2007l : iGrossi &: Springell l2009l ) . More mas- 
sive clusters than the one reported in this paper may be found in shallower X-ray observations down 



to 10 ^'^ergs s ^ cm ^, but one needs to survey an order of 5, 000 square degrees ([Finoguenov et al 



2009bl ) . which may be difficult to follow-up photometrically. However, a similar yield in the number 



of clusters can be achieved by reaching 10~ ergs s~ cm~ over 100 square degrees. Therefore, 
deep X-ray surveys at high spatial resolution have a unique window for cosmological studies at the 
time when Universe was only a third of its present age. Our observation also shows that detailed 
analysis of galaxy populations at this redshift regime is feasible with the current near-IR facilities. 
This is an encouraging result and motivates us to push our X-ray cluster survey forward. 

This study is based on data taken at the Subaru Telescope operated by the National Astronom- 
ical Observatory of Japan through program S09B-042. We thank Ichi Tanaka for his support during 
the observation. We thank the anonymous referees for their comments, which helped improve the 
paper. 
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Fig. 1. — Sample spectra of z = 1.6 galaxies. The dashed lines are the noise spectra. Some of the 
most prominent spectral features are indicated with the dotted lines. Zconf = and 1 mean secure 
and possible redshifts, respectively. The bottom spectrum is smoothed over 5 pixels. 
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Fig. 2. — i?z3.6/xm pseudo-color image of the twins. The image is 5 arcmin on a side. The contours 
show diffuse X-ray emission at the 1.5a, Scr, and 4.5o" levels in the 0.3-2.0 keV band and point 
sources are removed. The arrows indicate photo-z selected potential members (1.4 < Zphot < 1-8). 
The circles s how our objects with spectroscopic redshifts and the triangles indicate spectroscopic 
objects from iPapovich et al.l (|20ld ). The redshifts with '?' are possible redshifts. 
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Fig. 3. — z — K plotted against K. Galaxies in the left panel are selected first by photometric 
redshifts (1.4 < Zphot < 1.8) and then the remaining fore-/background contamination within the 
photo-2: slice is statistically subtracted. The dot and points show subtracted and remaining galaxies, 
respectively. We also show spectroscopic objects at 1.62 < Zgpec < 1-65 and the meanings of the 
symbols are shown in the panel. The dashed lines show the K-hand magnitude cut and 5cj limiting 
colors. The dotted lines are the model red sequences formed at Zf = 2, 3, 4, and 5 from bottom to 
top. The solid line is a biweight fit to bright red galaxies {K < 22 and z — K > 2.5). We apply 
only the statistical field subtraction without using photo-z in the right panel. 
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Fig. 4. — LFs of red galaxies. The top panel show the Ngai (after the field subtraction) and 
the bottom panel show Ngai normalized at ~ m*^ at each redshift. The top axis shows i^T-band 
magnitudes for the z = 1.6 galaxies, while the bottom axis shows -ftT-band magnitudes relative to 
at each redshift. The meanings of the symbols are shown in the plot. At ^ = 1.6, is derived 
by evolving m*^ at z = back in time assuming Zf = A. The points are sh i fted h orizontally to 
avoid overlapping. The LFs at z = and z = 1.1 are taken from lTanaka et aP ( 2008 ). The dashed 
line shows the magnitude cut of -R' = 23 for the z = 1.6 sample. 




Fig. 5. — The observed spectrum of the brightest red galaxy with secure redshift (dark and hght 
coded are the spectrum and its error) overlaid with the best-fitting model spectrum. The shaded 
region is not used in the fit to avoid effects of strong atmospheric extinction. 



